Behavior of a Magnetic Impurity in Graphene in the Presence of a Vacancy 
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With quantum Monte Carlo methods, we investigate the consequences of placing a magnetic 
adatom adjacent to a vacancy in a graphene sheet. We find that instead of the adatom properties 
depending on the energy of the adatom orbital, as in a single impurity problem, they develop a 
dependence on the energy of the split localized state associated with the single vacancy problem. 
Shifting the chemical potential through this experimentally more accessible energy scale reveals 
novel behavior in the spectral density, magnetic susceptibility, and the correlations of the adatom 
spin and charge with those of the conduction electrons. In general, the behavior of the adatom in 
the presence of a vacancy differs significantly from its behavior in the absence of a vacancy. 
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I. INTRODUCTION 

The interest in the consequences of defects in graphene 
has been mushroomingii"— Graphene's linear density of 
states imparts unusual properties to the behavior of its 
defects which in turn generates unusual physical proper- 
ties for the graphene sheets. The interest is especially 
high in controlling these defect-induced properties, par- 
ticularly for those cases where the adatom develops a 
magnetic moment. 

Usually vacancies or magnetic adatoms in graphene are 
individually studied, and they have been found to exhibit 
unusual properties. For example, a vacancy leads to the 
formation of a sharp doubly degenerate localized state 
sitting at the junction of the upper and lower bands. This 
localized state is quite robust to the addition of many 
vacancies. In a nearest-neighbor, tight-binging particle- 
hole symmetric model of graphene, the amplitude of the 
resonant state wavefunction is zero on all lattices sites 
that are on the same sublattice of the vacancy^ 

A magnetic adatom on graphene behaves quite differ- 
ently from the classic Anderson impurity physics coupling 
to a constant density of states. Generally, with a linear 
density of states, a Kondo effect docs not occur;^ and 
such observables as magnetic susceptibilities and spec- 
tral densities differ markedly from mean-field and con- 
stant density of state renormalization group theories. 
As a function of model parameters, the adatom mag- 
netic moment in graphene can cross-over from having 
free moment Curie-like temperature dependence to hav- 
ing a paramagnetic-like temperature independent one. 
The spectral densities lack a central peak. Instead, they 
have two peaks, one in the upper band and one in the 
lower, separated by roughly U, the Coulomb interaction 
associated with the impurity orbitali^ii 

In this work, we report the consequences of an adjacent 
vacancy-adatom pair. We find that when paired these de- 
fects exhibit new unusual features that differ from either 
the single-vacancy or the singie-adatom case. For exam- 
ple, the robust central resonance of the single vacancy 
problem disappears, even in the particle-hole symmetric 




FIG. 1: (color online), (a) The vacancy-adatom system. The 
vacancy is at the site Rua and the adatom (red ball) is at its 
nearest-neighbor site Rit. The next-nearest-neighbor site is 

Rla. 



case. The upper and lower peaks of the single adatom 
problem are each split with the two outer peaks of the 
now four separated by roughly U, and the positions of 
the inner two have only a weak dependence on U . These 
peaks are the doubly degenerate localized vacancy states 
shifted to the upper and lower bands by the presence 
of the adatom. The significance of this split is that by 
changing the chemical potential of this two-defect sys- 
tem, we can turn the adatom magnetic moment on and 
off by varying it over a energy range set by the localized 
vacancy state instead of a range set by the location of 
the adatom impurity orbital. We also find that the two- 
defect wavefunction has a small amplitude on sublattice 
sites the same as that of the vacancy, just as in the case 
of a single vacancy. 



II. MODEL AND METHODS 

Our starting point is the Anderson impurity model 
which has a single impurity orbital of energy ej, and 
Coulomb repulsion U inhibiting its occupancy by two 
electrons)^ The impurity orbital is coupled to a free- 
electron conduction band with hybridization of strength 
V. The total Hamiltonian is 

H = Hq + Hi + H2. 
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As standard, Hq is a tight-binding Hamiltonian which for 
graphene is 

Ho = -t [al^bj^+b]^ai^]-fi ^ [aj^ai^+bj^b.^] , 

<ij>,a <ij>,a 

where and bl^ create an electron with spin a at sites 
Ria and Rib on the A and B sub-lattices of graphene's 
hexagonal structure. In graphene the hopping matrix el- 
ement t > and is about 2.8 eV— . ^ is chemical potential 
and in ideal graphene its value is 0. Hi is the impurity 
Hamiltonian 

where d\ creates an electron with spin cr at the impu- 
rity orbital. Finally, H2 describes the hybridization be- 
tween the adatom impurity and one of graphene's carbon 
atoms. As shown in Fig. [TJ the tight-binding portion of 
the Hamiltonian is perturbed by creating a vacancy at 
site i?oa- In what follows wc place the adatom either on 
the opposite sublattice at or on the same sublatticc 
at Ria- For an adatom on site 

Our principal computational tool is the single- impurity 
quantum Monte Carlo (Hirsch-Fyc) algorithm for com- 
puting the thermodynamic properties of the adatom, and 
the method of Baycsian statistical inference for com- 
puting its spectral density. The Hirsch-Fye algorithmic 
naturally returns the imaginary-time Green's function 
Gd{T) = ^^Gda (t) of the impurity. With this Green's 
function, we can easily compute for the adatom orbital 
such basic quantities as the expected values of the total 
charge Ud, the square of the spin to^, and the doubly 
occupancy rid-^ndi- We note that 

ml^Ud- 2nd^ndi (1) 

A non-zero value of to^ indicates the formation of a mo- 
ment on the adatom orbital. The closer this value is to 
one the more fully developed is the moment. Wc also 
calculate the static adatom spin susceptibility 

X= / dT{md{T)mdm, (2) 

where P = T'^, mdir) = e^"mdiO)e-^" . 

Computing the imaginary-time Green's function also 
enables us to compute the spectral density = 
A(j(cj) by numerically solvingii 

^ ' J e-P'^ + 1 
—00 
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FIG. 2: (color online), (a) as a function of ^. (b) 
ndupnddovin as a function of /i. (c) as a function /i. (d) x 
as a function of /i. In all plates, U = O.SOf, Sd = —U/2, and 



The detailed procedure for doing this is presented in 
Ref. [ill. This Bayesian inference procedure is colloqui- 
ally, but imprecisely, also called the maximum entropy 
method. 

Invoking the bipartite nature of the lattice , which 
remains true even with the vacancy and adatom, 
and using the standard particle-hole transformation 
on one-sublatticc. we can prove that when Ed = 
—U 12. Gd{T,fi) = Gd{—T,—^) and consequently that 
Ao.(/i,a;) = A„{^^^—uj). These two results in turn im- 
ply symmetries with respect the sign of ^ in the various 
thermodynamic quantities of interest. These symmetries 
also mean that without loss of generality we can restrict 
our attention to the behavior of the system when ^ < 0. 

Wc calibrated our quantum Monte Carlo results by 
performing exact diagonalization studies^ of an 11-site 
graphene structure with a vacancy and adatom. This sys- 
tem is half-filled when there arc 12 electrons. We started 
with this filling and then systematically removed elec- 
trons, keeping the number of up and down electrons equal 
when their total number N was even or having one more 
up than down if their total number was odd. For the same 
system, we also performed quantum Monte Carlo calcu- 
lations and found excellent agreement in the computed 
thermodynamic quantities and correlation functions. 

Wc comment that the exact diagonalization calcula- 
tion is at zero temperature and in the canonical ensem- 
ble. The quantum Monte Carlo simulation is at a finite 
temperature and in the grand canonical ensemble. To 
compare quantum Monte Carlo with exact diagonaliza- 
tion, we did the simulation at a very low temperature 
(/3 = 52) and adjusted the chemical potential so that the 
average number of electrons closely matched that of the 
zero temperature systemJ^ 
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FIG. 3: (color online), (a) The spin susceptibility x with 
negative /.i as a function of T^^ (b) x with positive jj.. U — 
0.804, V = 0.50t and = -0.40t. 
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III. RESULTS 



A. Basic Thermodynamic Quantities 



FIG. 4: (color online). The spectral density A{iu) versus uj. 
(a)-(c) V = 0.504 and (d)-(f) V = l.Ot. 



In Fig. [2] we present our study of the behavior of sev- 
eral basic thermodynamic properties as a function of fi 
when the adatom is at Rib- Here = —U/2, so when 
^ = 0, the system possesses particle-hole symmetry. This 
symmetry fixes the total charge on the adatom at one. 
As we see from Fig. [2^, by shifting fj, lower, we reduce 
this charge. Figure [It shows the expected value of the 
adatom's double occupancy. The double occupancy is 
relatively small and becomes smaller as the value of /i is 
lowered, but increases as V increases. From these two 
curves we could predict the values of the rn^ from Eq. [1] 
We actually compute it independently. We show this 
measure of a local moment in Fig. Despite the gen- 
eral reduction of the total charge and double occupancy 
as we lower n, we see that for V = 0.50t the moment 
first increases and then drops in value. For V = 0.75t the 
moments arc much smaller with an even smaller initial 
increase. When V ~ l.Ot, the moment is even smaller, 
and we sec no initial enhancement. Fig. [2}1 shows the 
adatom's spin susceptibility in Eq. In it is mirrored 
the behavior of rn^. The enhancement for V = 0.50< is 
significantly larger than that for the other two cases. 

The principal differences among the results in Fig. [5] 
are for the V = 0.50t case initially decreasing more 
rapidly and varying differently as ii varies over the inter- 
val (-0.10f,-0.20i). We note that around ^ = -0.20t, na 
for V = O.SOi even becomes larger than that of = 0.75t. 
These results point to an anomalous situation, where as 
we decrease the total charge on the adatom orbital, which 
must be one or less, we increase the instantaneous net 
spin of these electrons. 

We comment that if were to place the adatom at Ria, 
we would find much smaller enhancements if any occur. 
We recall that in the absence of an adatom the wave- 
function of the localized vacancy state is non-zero at i?if, 
but zero at Ria- We therefore expect differences in the 
adatom properties for the two different placements. 

Figure |3] shows the adatom's spin susceptibility x ^-s 
a function of inverse temperature for different values of 
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FIG. 5: (color online). The spectral density A{ij) versus uj 
with different values of U. Here V = l.Ot and Sd = —U/2. 



positive and negative for V = 0.50t and Ed = —U/2. 
To numerical accuracy, the results for the different signs 
of /X are identical. For = 0, a particle-hole symmetric 
case, we see the susceptibility becoming temperature in- 
dependent at low temperature, even though a fairly well 
developed moment exists (Fig. [5J:). This co-occurrence 
is precedented but unusual. The temperature indepen- 
dence develops because the vacancy and adatom states 
are well below fi. As we lower /i a Curie-like tempera- 
ture dependence develops and then crosses over to a very 
weak temperature dependence. 



B. Spectral Densities 

We also calculated the adatom's spectral densities. For 
the V = 0.504 case we show our results in Fig. \^-c for 
^ = -0.054, -0.104 and -0.214. The dashed lines mark 
the location of the chemical potential: ui = ^. In all 
sub-figures U = 0.804, Ed = -U/2 = -0.404, T = 324"^ 
We clearly sec the four-peak structure noted earlier plus 



4 





/ 0.000 

/ (a) 

11 -0,002 
il -0.004 

— A— i.-.n 5m 






' ^ ' 0123456789 







/ 0.01 0. — 
IB 0.005. 

Jl (b) 

// -0.005. 






' 0123456789 10 



01234 5j 6789 10 

FIG. 6: (color online). (a) Charge-charge correlation function 
Ci (b) Spin-spin correlation function Si. U = O.SOi, V = 0.50t 
ea = -OAOt and T"^ = 32t"\ The insets show the tails of 
Ci and Si. 



the breakdown of particle-hole symmetry as soon as fi 
is slightly shifted from zero. Figure |4]i-f for V — l.Ot 
complements the V = 0.50i case. In both cases, as we 
lower fj, the former localized vacancy state peaks are first 
pushed away from w = and then are drawn back. Con- 
sequently, the dynamics of the local moment is initially 
controlled by the localized vacancy state and then by the 
adatom impurity state. 

We note that when /x is far from zero, the larger V 
results give a wider width to A{uj) peaks associated with 
the impurity state. This is not unexpected, but we also 
note that our procedure for determining the spectral den- 
sity tends to broaden high frequency fcatures^ii 

In Fig. [5l we study the spectral densities varying U 
with the hole-particle symmetry, Sd = —U/2 and /i = 0, 
and here T = 32t~^. In bulk transition metal, the value 
of U is about 5eV to lOeVr^ and in the case of transition 
metal atom in carbon-based materials, the value of U is 
about 2eV to 5eV.— We change the Coulomb interac- 
tion strength in the range (0.8t,3t), which overlaps the 
above ranges. We see that when U is increased, the outer 
two peaks become farther from each other while the two 
inner become closer, and we also note that these changes 
are small in this range. We expect that the two inner 
peaks will continue to move toward each other when U 
is increased. 

Wc comment that the strong central peak characteris- 
tic of a vacancy in graphene is absent when the adatom 
is present. In fact it becomes absent as soon as particle- 
hole symmetry is absent. The inner peaks in Fig. 2] are 
these fi = localized states shifted by the symmetry 
breaking and the presence of the adatom. With exact 
diagonalization we computed the local density of states 
for several such cases in the absence of the adatom. The 
only presence of the adatom perturbs the peak locations 
slightly. 



C. Correlation Functions 

The expected physics of the single impurity Anderson 
model in the absence of a vacancy and the Kondo effect 
is a double peak structure in A{uj) where the lower fre- 
quency peak is associated with a singly occupied state 
and the upper peak with an excitation to a doubly oc- 
cupied state. The single electron in this impurity state 
is anti-ferromagnetically correlated with the conduction 
band electrons and fluctuates in its spin polarization so 
that {nd-\—ndi) = 0. Fluctuations in charge are relatively 
weaker. 

An extension of the Hirsh-Fye algorithniii makes the 
computation of the the spin-spin Si and charge-charge 
Ci correlations between the adatom and conduction elec- 
trons possible. We show typical results in Fig. [51 Co and 
Sq are the on-site correlation functions for an adatom at 
Rib. The locations with even index are positions on the 
same sublattice as Rn, and those with an odd index are 
positions on the opposite sublattice. We see that both 
the spin-spin Si and charge-charge correlations Ci are 
relatively short-ranged. 

The plots of Si show that the adatom spin is always fcr- 
romagnetically coupled with the conduction electron spin 
at neighboring sites. For the single adatom case, this cor- 
relation shifts from ferromagnetic to anti-ferromagnetic 
as we shift /i below zero. In addition, the correlations 
are longer-ranged and less damped. We also note that Si 
as a function of displacement from the vacancy quickly 
becomes anti-ferromagnetic when the adatom spin and 
conduction electron spins are on the same sublattice but 
represents uncorrelatcd spins if they are on opposite sub- 
lattices. Presumably, as for the single vacancy case, the 
two-defect wavefunction is very small at the sites in the 
second case and larger at the sites in the first case. We 
further note that as we lower /x the spin correlations be- 
come monotonically suppressed, even though over the 
range of the shift where we are increasing and then de- 
creasing the adatom moment (Fig. This suppression 
is the opposite of the enhancement we observed for the 
single adatom case<i 

The charge-charge correlation function Ci is larger if 
the adatom and conduction electron charges are on the 
same sublattice. When we shift /i below zero, its ampli- 
tude on these sites also decreases. When we shift fi, the 
occupancy of the impurity orbital and conduction band 
shift from half filling so their charge exchange (the pos- 
sibility of hopping) enhances, and the expected value of 
the charge decreases (Fig. which in turn decreases the 
size of the correlations. This behavior is also opposite of 
what we observed for the single adatom caseii 

IV. REMARKS AND CONCLUSIONS 

Wc expect the strong features in A(aj) to appear in 
the system's density of states. Accordingly, the position 
of the chemical potential relative to the position of the 
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peaks in A(lo) can influence the physics of the adatom. 
When the chemical potential is near or at the location of 
the shifted localized vacancy state, the physics is that of 
a localized state. When away, the physics changes. Our 
results are consistent with this expectation. 

Our findings are a result of the impurity orbital now 
sharing its electrons between two states with two differ- 
ent energies instead of its electrons being associated with 
a single energy state. Still, these states must share a total 
charge of one or less, depending on the chemical poten- 
tial, and both must respect the energy penalty for dou- 
ble occupancy. In turn, these two states hybridize with 
each other through their hybridization with the conduc- 
tion electrons and exchange electrons. As we lower the 
chemical potential and begin to deplete the electron oc- 
cupancy in the upper localized state, at the same time 
reducing the double occupancy (Fig. [2]), we can increase 
the fluctuations of net non-zero spin and hence increase 
771^. When we leave the localized state behind, we are left 
with just the impurity state, which now has fewer remain- 
ing electrons and strongly suppressed spin and charge 
fluctuations. In general, as a function of the chemical 
potential, the presence of the vacancy causes the adatom 
impurity to behave differently, and in the case of its spin 
and charge correlations with the conduction electrons, 
for example, the behavior trends as a function of ^ are 
opposite to those of the single adatom case. 

It is clear that by shifting the chemical potential 
away from the half-filled particle-hole symmetric case 
we can have a magnet moment that we can enhance 



or diminish by changing the chemical potential. We re- 
mark that vacancies in graphene are readily generated 
by irradiation . We also note that organic groups ab- 
sorbed on the graphene sheet can generate the localized 
states similar to those induced by vacancies)^ So the 
organic groups also can be regarded as a kind of ideal va- 
cancies in graphene. And finally we remark that scanning 
tunneling microscopy (STM) opens the possibility for ex- 
perimentalists to place a variety of magnetic adatoms at 
precise locations on a graphene sheet.— In principle, our 
predictions are testable experimentally. 

We have not studied how our predictions change as the 
separation between the two defects increases, other than 
to show that the effects are most pronounced if the va- 
cancy and adatom are on different sublattice positions. 
Clearly, if the system has many of either or both defects, 
the results will likely change significantly. Many vacan- 
cies, for example, fill in graphene's pseudo-gap, removing, 
on the average, the linear density of states. 
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